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Abstract 

By manipulating the clustering coefficient of a network without changing its degree distribution, 
we examine the effect of clustering on the synchronization of phase oscillators on networks with 
Poisson and scale-free degree distributions. For both types of network, increased clustering hinders 
global synchronization as the network splits into dynamical clusters that oscillate at different 
frequencies. Surprisingly, in scale-free networks, clustering promotes the synchronization of the 
most connected nodes (hubs) even though it inhibits global synchronization. As a result, they show 
an additional, advanced transition instead of a single synchronization threshold. This cluster- 
enhanced synchronization of hubs may be relevant to the brain which is scale-free and highly 
clustered. 

PACS numbers: 05.45.Xt, 89.75.Hc 
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A keyproblem in the study of networks is the relation between network structure and 
function QQ]. Among the attributes most frequently used to characterize network structure 
are degree distribution and clustering coefficient. The latter, defined as the number of trian- 
gles divided by the number of connected triples , quantifies the tendency of neighbors of 
a node to be also neighbors of each other. We refer to this property as structural clustering 
to distinguish it horn, frequency or dynamical clustering jslQljal. 

A recent algorithm p allows one to manipulate the clustering coefficient by re- wiring the 
network, without changing its degree distribution. Kim used this technique to study the 
effect of clustering on the performance of variously structured Hopfield networks and we 
use it here to study its effect on synchronization. 

In search of factors that control the synchronization of networks of oscillators, re- 



searchers have studied the effects of different structural and statistical attributes j 
These studies have included systems of coupled mapsjoj ^jjsl 3|, continuous-time chaotic 
oscillators jl^ jl^ . spiking neurons jig, and phase oscillators [3] lfij \v\ fl^. Unsur- 
prisingly, shortcuts in small-world networks tend to improve synchronization compared to 
regular lattices 3]- However, other factors such as degree heterogeneity . 
maximum betweenness centralityQ], asymmetry and weighting of couplings 20 ] j^lj also 
play a role. 

In the present work we used Kim's procedure [5] to study the effects of structural clustering 
on the synchronization of networks of Kuramoto-like [2^ phase oscillators with Poisson and 
scale-free degree distributions. For both network types we found that increased clustering 
impedes global synchronization and magnifies the fluctuations of the global order parameter. 
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Clustering also qualitatively changes the onset of synchronization. At low clustering 
the transition is similar to that in the mean-field or globally coupled Kuramoto model, where 
a single subset of oscillators becomes entrained at a central frequency. With increasing 
coupling, this synchronized subset entrains larger portions of the remaining oscillators. At 
higher clustering, on the other hand, synchronization begins with more than one subset, 
each synchronized at a different frequency, and continues at increased coupling through the 
recruitment of remaining oscillators by the different subsets. We refer to these subsets as 
frequency clusters. In contrast to Poisson networks, characterized by a single synchronization 
transition, we find that strongly clustered scale-free networks show a second synchronization 
transition which is advanced relative to the principal transition, i.e. it occurs at a lower 
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value of the coupling strength. Thus, structural clustering promotes synchronization at low 
coupling, while inhibiting it at higher coupling. The advanced synchronization begins with 
the hubs, or highest- degree nodes. These results appear to be relevant to natural networks, 
such as the braiu H and g aa g .ia Q Q Q that have higher clustering eoeffloeuts thae those 
predicted by simple growth models. In the brain the distribution of functional connections 
and the probability of finding a link vs. distance are both scale free 1231. 

finnn n 

In contrast to approaches 11] 12] [14] [19] based on the master stability function |25(, we 
consider non-identical oscillators and their full dynamics, including states where only some 
of the oscillators are synchronized. 

We consider networks of oscillators obeying the coupled differential equations 

~dt = Ui + jk^ aij sin ^ ~ ^ 
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where < fa < 2tt are N phase variables, cjj are the randomly and uniformly distributed 
intrinsic frequencies, A is the coupling strength, and is the adjacency matrix (1 if i and 
j are connected, otherwise). The coupling strength is scaled by the average degree (k) of 
the nodes, averaged over the whole network. In the globally coupled Kuramoto model [22] 
where a^- = 1 , Vz ^ j, this reduces to scaling by the number of oscillators N. The above form 
makes the coupling strengths symmetric and weights all links equally[17], but in the case of a 
nonhomogeneous degree distribution, some oscillators may receive a stronger synchronizing 
signal because they have more neighbors. 

Results are for networks with N = 1000 nodes and average degree (k) 20. Quali- 
tatively similar results were obtained with N = 5000 and (k) = 6. We consider a random 
(Poisson) network and a scale-free network generated by the Barabasi-Albert preferential 
attachment algorithm [3], and families of networks derived from each of these by changing 
the clustering coefficient 7. To vary 7, we use a stochastic rewiring algorithm^ that rear- 
ranges connections with a bias toward increased clustering. The procedure is as follows: 1) 
Pick randomly two existing links. 2) Compute whether interchanging these links increases 
or decreases the total number of triangles. Perform the interchange only if it increases that 
number. 3) Repeat these steps until the desired clustering is achieved. (One can reduce 
the clustering coefficient by reversing the acceptance criterion). Since this algorithm only 
rewires connections and does not change the degree of any node, the degree distribution as 
well as the degree sequence is fixed. 



The random uoi values were uniform distributed over the interval 0.9 < Ui < 1.1 and the 
initial phases were also random. We ran the dynamics at a series of increasing values of 
A, integrating the equations using a simple Euler method with step size 0.02. After 100 
time units to allow for relaxation to a steady state, we measured the time-averaged (over 



500 time units) synchronization order parameter! 

(2) 



m = 



I 3=1 
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m is of order 1/y/N if the oscillators are uncorrelated and approaches 1 when all are in 
phase. In addition, we measured the frequency Q{ = j^j jof each oscillator over the same 
500-unit time interval after relaxation. The brackets {} signify time averaging. This 
measurement reveals the collective behavior in detail. 

To study the effect of clustering on the transition to synchronization, we first examine 
the order parameter m as a function of the coupling strength A. Figure ^ shows plots 
of m vs. A for networks with Poisson and scale-free degree distributions at several values 
of the clustering coefficient 7. In each plot the data were averaged over time, over several 
realizations of the intrinsic frequencies and over several network rearrangements. The degree 
sequence, however, was the same in all cases. Poisson networks (fig.lA) behave more simply 
than scale-free ones. The natural (low-clustering) random network shows a rapid transition 
to order above A = 0.12. Increasing clustering inhibits global synchronization and makes 
this transition less steep. The lowered mean values of m in the highly clustered networks 
are associated with extremely large temporal fluctuations around the mean. 

In scale-free networks (fig. IB, C), clustering has a more complicated effect. The network 
with 7 = shows the onset of synchronization (upward turn in the graph) near A ~ 0.07 and 
the transition to full synchronization is not as steep as in the Poisson case. Surprisingly, 
increased clustering slightly enhances the order parameter at weak coupling but suppresses it 
at strong coupling. For intermediate 7 this results in two separate transitions, an advanced 
transition (threshold at A ~ 0.03) to partial synchronization followed by a delayed one 
(A ~ 0.15) to full synchronization. 

For more insight into the synchronization transitions we examine scatter plots of observed 
frequency fi, vs. intrinsic frequency cjj. Each plot in figures EHU shows the population's 
behavior for one realization of the random variables. In the absence of interactions, all 
points would fall along the line Qi = o^. In the globally coupled (mean-field) case[22j, the 




FIG. 1: (color online) Order parameter m vs. coupling strength A, for different values of the 
clustering coefficient 7. (A) Poisson degree distribution. (B,C) Scale-free degree distribution. 
(C) Close-up of the transition region showing that increased clustering leads to an advanced (lower- 
A) transition. 



nodes differ only by their intrinsic fequencies. Therefore, the observed frequency is in every 
case a single- valued function Q(u) of the intrinsic frequency, whose shape changes with A. In 
particular, the curves begin to flatten around uo ~ 1, as oscillators with frequencies closest to 
the average are the first to synchronize. The flat portion broadens with increasing coupling. 
For Poisson networks (fig. EI) at low clustering, the same qualitative picture holds as in the 
mean-field case except that there is more scatter and the curve is blurred. This is due to the 
fact that each oscillator's behavior now depends not only on its intrinsic frequency and the 
global average, but also on the details of its local neighborhood. As in the mean-field case, 
the oscillators near the extremes of the u distribution are the last to become fully entrained. 
Unlike the mean-field case, some remain unentrained even when the synchronized group is 
well-established and has recruited members from the entire intrinsic frequency range. At 
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FIG. 2: Poisson networks: Scatter plots of average oscillation frequency f2 (vertical axis) vs. 
intrinsic frequency oo (horizontal). At zero coupling, the points all lie on the line Q = u. With 
increasing coupling, more oscillators line up at synchronized frequencies. For high clustering, 
several bands form at different frequencies, rather than one frequency as in the low-7 case. 



higher clustering, the scatter plots differ even more from the mean-field case, and they do 
so in two ways. First, the scatter is more pronounced. Second, the horizontal striations 
indicate that the oscillators cluster in subgroups (frequency clusters) oscillating at different 
frequencies. As A increases, the frequency clusters converge and eventually merge, but full 
synchronization requires a stronger coupling compared to networks with low clustering. The 
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FIG. 3: Scale-free networks: Scatter plots of O vs. lo. Synchronized subsets (horizontal bands) 
begin to form at weaker coupling for the more clustered networks. Frequency clustering is evident 
at high 7. 



beating of the different frequencies accounts for the low value of the global order parameter 
and its large fluctuations. 

In the scale-free case, the scatter plots (figure EJ) confirm that the highly clustered 
networks begin to synchronize at a weaker coupling than the less clustered ones. Initially, 
however, this synchronization only affects a subset of the nodes, while the remaining nodes 
still fall close to the line Q = uj as if they were not interacting at all. This is in contrast 
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FIG. 4: Scale-free networks: Average frequency £1 vs degree at several values of clustering 7 and 
coupling strength A. The degree (horizontal axis) is plotted on a logarithmic scale. Synchronization 
begins with the hubs (nodes with the highest k) and progresses downward to those with lower k. 



to the Poisson scatter plots which begin to flatten near the center of the frequency range 
as the transition begins. Since scale-free networks have a strongly heterogeneous degree 
distribution, the entrainment of an oscillator depends strongly on the number of its inputs. 
When Q is plotted against the degree k of each node as in figure El it is apparent that 
the higher-degree nodes (hubs) begin to synchronize first, while the lower-degree nodes 
synchronize as the coupling continues to increase. One similarity with the Poisson case 
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is that structural clustering enhances frequency clustering. While clustering promotes the 
formation of synchronized frequency clusters among the hubs, it inhibits the synchronization 
of the network as a whole. The early hub synchronization accounts for the slightly enhanced 
order parameter at weak coupling. 

In conclusion, we examined the synchronization of networks of non-identical coupled 
phase oscillators with both Poisson and scale-free degree distributions, and studied the 
effects of varying the clustering coefficient without affecting the degree distribution. Our 
first main result, for both types of networks, is that clustering encourages the formation 
of sub-populations synchronized at different frequencies (frequency clusters), and thus dis- 
courages full global synchronization at a single frequency. Our second key finding concerns 
the scale-free case, where we found that, despite the increased difficulty of full synchroniza- 
tion, higher clustering actually promotes the onset of partial synchronization of the hubs. 
Scale-free networks with high clustering thus appear to undergo two separate transitions as 
the coupling strength increases: an early transition to partial synchronization and a delayed 
one to full synchronization. The first transition involves only the hubs while leaving the 
majority almost unaffected. The hubs seem to form the growth nuclei for the synchronized 
state. Interestingly, they are able to synchronize even though many of their inputs come 
from lower degree nodes which are not yet synchronized. It is known that in ordinary, 
low-clustered scale-free networks the relaxation time for synchronization of hubs is shorter 
than that of less connected nodes, but this does not account for the effect of structural clus- 
tering in promoting their synchronization. It is possible that increasing clustering may, as a 



by-product, increase the assortativity of degree mixing [291]. so that hubs become more likely 
to connect to each other. In other models of clustered networks it was found that triangles 
preferential* .nc.ude ai g her- d e gre e node*,. Cluing may als0 affect the bet™ 
centrality distributions. 

This advanced transition may play a role in highly clustered networks, for instance in the 
human brain. In natural as opposed to abstract networks, the cost associated with long- 
range connections often gives them a tendency toward clustering. For the ganglion of C. 
elegans [24], the clustering coefficient is 7 = 0.28. For the human brain j^J, the connectivity 
is essentially scale-free in different regions, and clustering lies three orders of magnitude 



above that o 
in the brain 



equivalent random networks. This and the weak synaptic coupling strength 
3 [ill point to the potential importance of cluster-enhanced synchronization 
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of oscillating hubs. 
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